Testing the Copernican principle by constraining spatial homogeneity 
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We present a new program for placing constraints on radial inhomogeneity in a dark-energy dom- 
inated universe. For the first time, we quantify restrictions on violations of the Copernican principle 
without neglecting dark energy. This is important because any violation of this principle would 
interfere with our interpretation of any dark-energy evolution. In particular, we find that current 
observations place reasonably tight constraints on possible late-time violations of the Copernican 
principle: the Figure of Merit in the parameter space of amplitude and scale of a spherical inhomo- 
geneity around the observer has to be increased by a factor of 3 so as to confirm the Copernican 
principle. Then, by marginalising over possible radial inhomogeneity we provide the first constraints 
on the cosmological constant which are free of the homogeneity prior prevalent in cosmological mod- 
elling. 
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Introduction The Copernican principle states that 
humans are not privileged observers of the universe and 
provides our philosophical basis for assuming that on 
the largest scales the universe is spatially homogeneous. 
While it is one of the foundational aspects of modern 
cosmology, this assumption remains untested outside of 
the standard paradigm. Though it may seem pedantic to 
test something so obvious, the standard paradigm itself 
is built on shaky foundations, relying on an unexplained, 
gravitationally repulsive, dark-energy component for ob- 
servations to fit the model. The implications of this can- 
not be overstated. Assuming that the laws of physics do 
apply equally everywhere in the universe, the only non- 
copernican configuration possible is one in which we live 
in a place that originates from special initial conditions. 
The implications for our observations today would be an 
expansion rate that is sourced by different energy compo- 
nents (matter, curvature, etc.) in different directions and 
at different distances. As current observations of the ex- 
pansion rate of the universe suggest isotropy around us, 
any model that tries to explain part of the expansion his- 
tory as a result of spatial variations inevitably places the 
observer near the centre of a spherically symmetric mat- 
ter distribution. Because we can observe changes in the 
expansion of the universe as a function of redshift only, 
it remains hard to disentangle temporal evolution from 
possible spatial variations spherically around us. Indeed, 
there has even been active debate on models whereby 
(part of the) dark energy is replaced by radial inhomo- 
geneity about us, in violation of the Copernican principle. 
While this seems highly unlikely and challenging to make 
work (see [1 , 2] for reviews) , it highlights how intertwined 
are attempts to detect any evolution of dark energy to 
the Copernican principle. 

To rule out spherically symmetric models, consistency 



tests have been proposed (see [2] for a review), but all 
suffer as they lack at present a quantifiable measure of 
deviations from homogeneity. Attempts have been made 
to measure the homogeneity of the universe, albeit ignor- 
ing possible violation of the copernican principle [3, 4], 
We offer here a new approach, by directly constraining 
any radial change in the density, Hubble rate and cur- 
vature assuming the Copernican principle is false and 
that dark energy exists. We model the universe spher- 
ically symmetric, containing both matter and dark en- 
ergy in the form of a cosmological constant. We provide 
marginalised constraints on the amplitude and scale of 
any spatial inhomogeneity, over- and under-dense, using 
all relevant available observational data. 

We quantify the deviation from Copernicanism by 
comparing the constraints we find on the spherical inho- 
mogeneity to the theoretical constraints on the existence 
of the same inhomogeneity in a Copernican universe. If 
the universe is Copernican, the temperature fluctuations 
observed in the CMB come from the same distribution as 
do matter perturbations around us. Observing the CMB, 
one can calculate the probability of having a spherical in- 
homogeneity of certain dimensions around us. The ratio 
of allowed volumes in parameter space, on the one hand 
constrained by observations on the lightcone and on the 
other hand by the distribution of the distant CMB den- 
sity fluctuations, gives a measure of how far we are from 
establishing the Copernican principle observationally. 

Finally, by marginalising over the inhomogeneity pa- 
rameters, we provide constraints on A which are indepen- 
dent of the homogeneity assumption and any reliance on 
the Copernican principle, thus being the most robust to 
date. 

The model We model radial inhomogeneity assum- 
ing a spherically symmetric Lemaitre-Tolman-Bondi so- 
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lution including a cosmological constant A (ALTB, see 
e.g. [5-7]). The metric is given by (c = 1) 



ds 2 



-dt 2 



1 — k{r)r 



: dr 2 + al(t,r)r 2 dn 2 , (1) 



where the radial (ay) and angular (a±) scale factors are 
related by an = (fltj_r)'. A prime denotes partial deriva- 
tion with respect to the coordinate radius r. The cur- 
vature k — k(r) is not constant but is instead a free 
function. The FLRW limit is k — >• const., and aj_ = an. 
The two scale factors define two Hubble rates: 



H ± =H ± (t,r) 
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The analogue of the Friedmann equation in this space- 
time is then given by H\ = M/a\ — k/a\ + A/3, 
where M = M(r) is another free function of r related 
to the locally measured matter density 8irG p m (t, r) = 
(Mr 3 ) 1 1 a\\a\r 2 , which obeys the conservation equation 
p m + (2H± + H\\)p m = 0. We introduce dimensionless 
density parameters for the CDM and curvature, by anal- 
ogy with the FLRW models: 
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so that Sl m (r) + Hfe(r) + ^(r) = 1. Note that in the 
previous equation the gauge fixing a±(to,r) = 1 has been 
used. Moreover, depends on r because the present- 
day expansion rate H± is inhomogeneous. Using (3) the 
Friedmann equation takes on its familiar form: 
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Integrating the Friedmann equation from the time of the 
big bang tbb(?*) to some later time t yields the age of the 
universe at a given (t, r): 



t — tbb — 
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Hence there is a relation between the functions ibb, ^fe 
and f2 m . Therefore the ALTB model is specified by two 
free functional degrees of freedom, and we will use ^(r) 
and tbb(f). By demanding a homogeneous age of the 
universe we fix the bang function to zero, ibb(r) = 0- 
This ensures the absence of decaying modes in the matter 
density [8, 9], in agreement with the standard inflationary 
scenario. 

We parametrize the only left freedom with the curva- 
ture function with the monotonic profile 



k a (r) = h + (k c - k b ) P 3 
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where is the comoving radius of the spherical inhomo- 
geneity and 

!1 for < x < a 

l„ e -|^f( 1 -f^f) f or a < x < i . 
for x > 1 

Now k a (r) is C n everywhere. We choose n = 3, such that 
the metric is C 2 and the Riemann curvature is C°. For 
r > rb the curvature profile equals k b such that there the 
metric describes exactly the ACDM model. The central 
over- or under-density, determined by curvature k c , is 
automatically compensated by a surrounding under- or 
over-dense shell. The smoothness of the transition from 
k c to fcf, is parameterized by a. Hence r^, k c , k b and a 
are free parameters. 

Finally, on the past light cone of a central observer, t(z) 
and r(z) are determined by the differential equations for 
radial null geodesies, 
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where Hu and ay are evaluated on the light cone. The 
area (cIa) and luminosity (d^) distances are given by 



d A (z) = a ± (t(z),r(z)) r(z), d L = (1 + z) 2 d 
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Data & Observables 

Hp: We compare a local average of the Hubble rate [10], 
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H±{r{z),t{z))dz, (9) 



with the value measured by [11] of -ffFrccdman — 72 ± 8 
kms -1 Mpc -1 . We choose z min = 0.005 and z max = 0.1. 
Supernovae la: We use the SNLS3 catalogue [12], 
which consists of 472 type la supernovae in the redshift 
range z = 0.01 — 1.39. 

CMB: We fit the CMB according to the method pre- 
sented in [13, 14], in which an effective FLRW metric 
is used to account for the different area distance to the 
surface of last scattering as compared to the homoge- 
neous background model. This method ignores isocurva- 
ture modes (consistent with the choice of a homogeneous 
bing bang) , assumes a standard number of relativistic de- 
grees of freedom and a standard power spectrum, all of 
which would change the constraints [2] . We fit our model 
to the WMAP 7-year data release [15]. 
BAO: The sound horizon during the drag epoch is im- 
printed in the galaxy correlation function. In an inho- 
mogeneous model this is an ellipsoid with proper scales, 
when viewed from the centre 
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where 6 S is the angle that the acoustic scale subtends 
on the sky, and z s is the redshift interval corresponding 
to the acoustic scale in the radial direction. The sound 
horizon d s is calculated assuming a homogeneous early 
universe. We use observations from SDSS, 6DFGS and 
WiggleZ [16-18], as compiled in [19]. 
Compton y-distortion: Off-centre observers see a 
large dipole in the CMB in an inhomogeneous universe. 
CMB photons are scattered from inside our past light- 
cone into our line-of-sight by off-centre reionized struc- 
tures which act as mirrors. The spectrum observed by 
the central observer is, therefore, a mixture of black- 
body spectra with different temperatures, producing a 
distorted black-body spectrum. In the single-scattering 
and linear approximations, and when the temperature 
anisotropy is dominated by the induced dipole /3, the y- 
distortion can be written as [13]: 
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where r is the comoving distance down the light cone 
and r le marks the reionization epoch. The time depen- 
dence of the optical depth is given by dr/dt = ot n e (t) = 
/b (1 — yke/2) p m (t)/m p , where ot is the Thomson 
cross section, /b = Pb/Pm is the baryon fraction, Yne is 
the helium mass fraction and m p is the proton mass. The 
dipole /3 is found by integrating the geodesic equations 
(7) in the negative and positive r-directions starting from 
an observer at {t(z),r(z)} back to the surface of last scat- 
tering. The difference in redshift between the two direc- 
tions is then approximately translated into the dipole ob- 
served by the scatter: /3 = zJ)j(2+z++zJ) . The 2a 
upper limit from the COBE satellite [20] is y < 1.5xl0 -5 . 
kSZ: The dipole /3 affects the observed CMB also 
through the kSZ effect [21]: hot electrons inside an 
overdensity distort the CMB spectrum through inverse 
Compton scattering, in which low energy CMB photons 
receive energy boosts during collisions with the high- 
energy electrons. We will focus here on the so-called 
"linear kSZ effect" [22], in which the effect due to all free 
electrons in the reionized universe is taken into account. 
Using the Limber approximation, the kSZ power at mul- 
tipole I is given by [23] 
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dimensionless 
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power spectrum of the background model and the func- 
tion k(r) = k(k(r),z(r)) is necessary to "isotropize" the 
angular and radial wave numbers which in an inhomo- 
geneous universe evolve differently: k(k,z) = k[(l + 
z)a x {t,r{z)f^ H {t,r{z))y s ]/[{l + z)a ± (z) 2 / 3 a|| (z) 1 / 3 ]. 
We constrain inhomogeneous models using a top hat 
prior < 1(1 + 1) [CgJaooo + Cfllooo] < 59/xK 2 , based 
on the results from SPT [24]. 




Figure 1. Marginalised constraints on rt and 5o from different 
combinations of data, C refers to CMB, F to Ho, B to BAO, S 
to SN, H to age data H(z), K to kSZ, Y to Compton-y distor- 
tion. The constraints are ordered from left to right and top to 
bottom, by the allowed area. Naturally all datasets combined 
give the strongest constraints. See Fig. 2 for a comparison of 
the CFBSHKY plot to the Copernican prior. 



Age data: Finally, we constrain radial inhomogeneity 
also by means of galaxy ages [25] using the maximum 
stellar age data from [26-28]. 

Copernican prior Given a gaussian density field, 
the root mean square of density perturbations inside a 
sphere of radius L around any point (hence also around 
the observer) today is 



T 2 

dk— 2 P(k) 



3ji(Xfc) V 
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where P(k) is the matter power spectrum today in- 
ferred from the CMB temperature spectrum, assuming 
a Copernican universe, as a function of wavenumber 
k, and ji is the spherical Bessel function of the first 
kind [29]. We calculate 07, for the radius L at which the 
central over/under-density makes the transition to the 
surrounding mass-compensating under/over-dense shell. 
This transition occurs at a radius L slightly smaller that 
TV Then we compare the actual density perturbation 
of a certain inhomogeneity So = M(r < L)/M — 1 and 

define the probability as P(Sq,L) cx exp —^(8q/(Tl) 2 . 

Note that we use the same normalization for different L 
so that the probability of having S — is the same. 

Results We present our results in terms of the non- 
Copernican parameters (*>o and rj,, marginalising over all 
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Figure 2. Marginalised posterior probability of r& and So from 
all data sets combined (blue contours) at 68%, 95% and 99% 
confidence level (c.L), compared to the Copernican prior ob- 
tained by fitting CMB only, (red to gray coloring) at 68%, 
95% and 99% c.L The ratio of the 99%-c.l. surfaces is roughly 
a factor of 3. Constraints on ri, and So need to improve by a 
factor of 3 to confirm the Copernican principle. 




0.65 0.7 0.75 



Figure 3. Marginalized posterior probability of the relative 
amount of dark energy (£2a) in the universe, when ignor- 
ing inhomogeneity (red dashed line), = 0.70 ± 0.04 at 
95% c.L, and when including inhomogeneity and marginaliz- 
ing over it (solid black line), Q. A = 0.69 ± 0.05 at 95% c.L 
Both constraints come from CMB, H , BAO, SN, H(z), kSZ 
and Compton-y. 



mogeneity which arises at late times, and do not apply 
to possible non-Copernican properties at early times. It 
is an important area for future study to generalise our 
results to these important cases. 



other aforementioned parameters. It should be clear that 
for a Copernican universe we have Sq = 0. In Fig. 1 we 
show the constraints on the non-Copernican parameters 
from a number of possible combinations of datasets, or- 
dered by constraining power. In all cases we use the 
CMB and Hq constrains in combination with at least 
one other dataset. In this combination it turns out the 
SNe are most constraining. Not surprisingly all datasets 
combined provide the strongest constraints, only allow- 
ing for a narrow range of density contrasts, however for 
many radii. Clearly when 8q — > 0, constraints on the 
radius allow for n, — > oo. 

The key result of this paper is shown in Fig. 2, where 
we compare the constraints from all datasets to the con- 
straints from purely the Copernican prior. We find that 
constraints on and 5q need to improve by a factor of 
3 to confirm the Copernican principle, because then ob- 
servations on the lightcone constrain inhomogeneity to 
within the range allowed for by the CMB power spec- 
trum. 

Finally in Fig. 3 we show for the first time constraints 
on the cosmological constant, A, marginalized over the ef- 
fect of inhomogeneities around us, compared to the same 
constraints without taking into account inhomogeneity. 
We find that error bars increase by 25% if one marginal- 
izes over inhomogeneity. 

Our constraints depend on a variety of assumptions, 
in particular that dark energy is described by a cosmo- 
logical constant [10, 30] and that general relativity is the 
correct theory of gravity. Our constraints apply to inho- 
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